A=1/eps;B=-eps;
C=[55,56,57,pi];
D=[pi,55,56,57];
S1=sprintf('%+15.5f',A)
S2=sprintf('%+.5e',B)
